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Abstract. We present Odin observations of the AGB star W Hya in the ground-state transition of ortho-ILO, lio - loi. at 
557 GHz. The line is clearly of circumstellar origin. Radiative transfer modelling of the water lines observed by Odin and ISO 
results in a mass-loss rate of (2.5+0.5) 10~ 7 M Q yr , and a circumstellar H 2 abundance of (2.0+1.0) 10~ 3 . The inferred mass- 
loss rate is consistent with that obtained from modelling the circumstellar CO radio line emission, and also with that obtained 
from the dust emission modelling combined with a dynamical model for the outflow. The very high water abundance, higher 
than the cosmic oxygen abundance, can be explained by invoking an injection of excess water from evaporating icy bodies in 
the system. The required extra mass of water is quite small, on the order of ~ 0.1 M e . 
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1. Introduction 

One of the main characteristics of an Asymptotic Giant Branch 
(AGB) star is its cool, low-velocity outflow. In this outflow, 
molecules are formed according to the relative chemical abun- 
dance of carbon and oxygen (Olofsson 2003 Millar 2003 and 
the references therein). In addition to H2, CO is an abundant 
species irrespective of the C/O ratio, but in C/O < 1 (O-rich) ob- 
jects its abundance is rivaled by that of H2O, while in C/O > 1 
(C-rich) objects HCN and C2H2 also have sizable abundances. 
Likewise, there are dust signatures of amorphous carbon and 
SiC when C/O is larger than unity, (Treffers & Cohen 119741 
while silicates are seen if the reverse is the case (Gilman l 19691 . 

Early AGB stars tend to be semiregular variables with a 
low mass-loss rate, and these evolve into long-period Mira-type 
variables where mass loss is more pronounced. The lifetime of 



a mass-losing AGB star is estimated to be about 10 5 yr, and be- 
fore evolving off the AGB the star undergoes an intense epoch 
of mass loss when most of the initial mass is lost, leaving be- 
hind a hot C/O core which then ionizes the surrounding ma- 
terial resulting in a planetary nebula (see e.g., Iben & Renzini 
IT3351 Habing [T996l . 

In this paper we discuss circumstellar H2O observations of 
the well-studied O-rich semiregular variable W Hya, with a pe- 
riod of 382 days (Lebzelter et al. 2000 1. It has been detected in 
a number of circumstellar molecular line emissions, e.g., CO 
(Wannier & Sahai 119861 Loup et al. 119931 Cernicharo et al. 
I1W71 Kerschbuam & 01ofsson[T333}, SiO maser (Pardo et al. 
2004 1 and OH maser lines (Szymczak et al. QP98 Chapman 
et al. IT994I Etoka et al. 120031 . It was also detected in water 
lines using the Infrared Space Observatory's (ISO) Short- and 
Long Wavelength Spectrometers (SWS and LWS) by Neufeld 
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Swedish Space Corporation was the industrial prime contractor and is 
operating Odin. 



et al. (1996 1 and Barlow et al. (1996 1, respectively. Recently, 
the Submillimeter Wave Astronomy Satellite (SWAS) detected 
the ground-state ortho-H20 transition at 557 GHz, 1 10 - 1 01, as 
reported by Harwit & Bergin ( 2002 1. 

O-rich semiregular variable stars are seen to have multi- 
ple dust species condensing in their outflows (Cami 2002 1. The 
10//m silicate feature, dominant in the spectra of Mira vari- 
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ables, may not be the strongest emission feature. This is clearly 
seen in the spectral energy distribution (SED) of W Hya. By fit- 
ting the SED with multiple dust species, Justtanont et al. (2004i 
derived a dust mass-loss rate of 3 10~ 10 M Q yr _1 (this is a lower 
limit as the 13/mi dust feature was not fitted due to the lack of 
a suitable identification). This suggests that W Hya is a low- 
mass-loss-rate AGB star. However, its dust envelope is known 
to be very large, extending out to 40' (Hawkins 1990 1. 



2. Observations 

The spectrum of the H2O ( 1 1() - loi) line at 557 GHz (A = 
538.29/mi) was observed by the Odin satellite towards W Hya 
using the autocorrelation spectrometer (Nordh et al. 120031 
Frisk et al. 2003 1. The resulting spectrum is a combination of 
two separate observing sessions, in December 2002 and July 
2003, of a total of 185 orbits. As the position-switching mode 
was used in these observations, the on-source integration time 
is 53.4 hours. The average system temperature for this fre- 
quency is 3 200 K. The calibration was done by switching be- 
tween the internal hot load and the sky. The frequency reso- 
lution is 1 MHz (= 0.54 km s _I at 557 GHz). A sinusoidal 
baseline was subtracted and the data rebinned. Fig.^shows the 
spectrum in terms of the observed antenna temperature, 7a, 
which has not been corrected for the beam efficiency of 90%. 

The peak antenna temperature is 0.05 + 0.02 K. The line 
flux is estimated to be (3.6+1.2) 10 -21 WcrrT 2 , after correcting 
for the aperture efficiency of 70% (Frisk et al. 2003 1. The ab- 
solute uncertainty in the line flux is large as the data are noisy, 
and close to the line there are features which are the artefacts 
caused by the ends and the joining of two correlator bands. 
These are seen as absorption features at Vlsr of 4, 23 and 30 
km s _1 (Fig[Q. The estimated uncertainty includes that due to 
the noise and the sinusoidal baseline subtraction, plus an abso- 
lute flux calibration uncertainty of 10%. The latter is obtained 
by comparing the line fluxes of Odin and SWAS data for a 
number of sources (Hjalmarson et al. I2003i . The line inten- 
sity is consistent with that obtained by SWAS after scaling the 
latter with the ratio of the Odin and SWAS beam sizes, 3.37 
(Hjalmarson et al. l2003> . 

The centre velocity and the expansion velocity are esti- 
mated to be 42 + 1.5 kms -1 (LSR scale) and 7.0+1.0 kms -1 , 
respectively. These values are consistent with the results from 
CO 7=1-0 data (Kerschbaum & Olofsson TF%% , 7=2-1 data 
(Cernicharo et al. lT997l . 7=3-2 and 4-3 data (Young [T995l . 



3. Modelling of the molecular line emission 

In this section, we present our models for both the H2O and 
CO line emission from the circumstellar envelope of W Hya. 
The models used are based on the assumptions of a constant, 
isotropic mass-loss rate, and a constant expansion velocity, 
which results in a spherically symmetric gas outflow with an 
r~ 2 density distribution. The mass-loss rate refers to the H2 
loss rate which is the quantity used throughout this paper. This 
value needs to be multiplied by a factor of 1 .4 to get the total 
gas mass-loss rate, due mainly to the presence of He. 
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Fig. 1. The ground-state transition of ortho-H20, 1 10 - lni, as 
detected by Odin towards W Hya. The three thick vertical lines 
mark where the two correlator bands overlap and end. 



In order to model the radiative transfer, we need to know the 
stellar luminosity. From the spectral energy distribution (SED), 
the peak of the spectrum corresponds to an effective tempera- 



ture of 2 500 K (Justtanont et al. 2004 ), which agrees well with 
that derived by Haniff et al. (I1995i . The Hipparcos distance, 
taken directly from the Hipparcos Catalogue (ESA 1997 1 par- 
allax, is 1 15 pc. The derived luminosity for this distance turns 
out to be very high - 1.16 10 4 L (estimated from a fit to the 
SED). However, the Hipparcos data were reanalysedby Knapp 
et al. (2003 1. The revised distance of 78 pc, results in a lumi- 
nosity of 5 400 L Q , more in line with it being an early- AGB 
star. The resulting radius is 2.73 10 13 cm. From this point on, 
we will assume a distance of 78 pc for our calculations. 

3.1. CO lines 

The code used to calculate the CO line emission is based on 
the work by Schoier & Olofsson ( 2001 ). It employs the Monte- 
Carlo method to determine the excitation in the lowest 40 ro- 
tational levels in the ground and first excited vibrational states. 
The collisional rates are taken from the CO-H2 rates by Flower 
(2001), and extrapolated to higher rotational levels and tem- 
peratures as described in Schoier et al. (2005 1 1 . The rates for 
ortho- and para-H2 collisions are weighted assuming an or- 
tho/para ratio of three. The outer radius is defined by the pho- 
todissociation of CO, and it is obtained following the results of 
Mamon et al. ( 1988 1. In order to convert the CO gas mass-loss 
rate to the total gas mass-loss rate, we assume a CO abundance, 
[CO]/[H 2 ],of210- 4 . 



1 The molecular datafile is available for download at 
www. strw. leidenuniv.nl/~moldata. 
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Fig. 2. The calculated best-fit model line profiles for CO 7=1- 
0, 2-1, 3-2, and 4-3 (continuous line) and the observed profile 
(histogram). 

In the fitting of observed line fluxes of CO 7=1-0 up to 
4-3, along with the ISO-LWS 7=16-15 and 17-16 lines, the 
gas kinetic temperature is calculated with a dust-drag heating 
balanced by the adiabatic and radiative cooling. The radiative 
cooling is provided by both CO and H2O molecules. The dust 
parameters entering in the heating term, due to collisions be- 
tween dust and gas particles, are taken to be consistent with 
those obtained from the dust modelling. In the case of a low- 
mass-loss-rate object such as W Hya, the CO molecules are 
mainly excited by the stellar radiation field and the temperature 
structure plays only a minor role (Schoier & Olofsson 2001). A 
best-fit model is obtained in terms of integrated line intensities 
and in the^ 2 sense, 
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where / is the total integrated line intensity, cr, is the uncer- 
tainty in observation i, and the summation is done over all inde- 
pendent observations N. The errors in the observed radio line 
intensities are assumed to be 20%, while the ISO lines have 
an uncertainty of 30%. The result is a stellar mass-loss rate of 
1.3 10 7 Moyr" 1 . 

Although our calculated line fluxes are close to the ob- 
served, the line profiles do not match. We consistently expect 
to see double-peaked profiles while the observations tell us 
otherwise. A way to remedy this is to decrease the envelope 
size and increase the input mass-loss rate. A mass-loss rate of 

1.3 10 _7 M o yr _1 results in a CO photodissociation radius of 

1.4 10 16 cm according to Mamon et al. fT9"88). However, our 
best-fit model, taking into account also the line shapes, has an 
outer CO radius of 8 10 15 cm with a mass-loss rate of 2 10" 7 
M yr _1 . The results of the best-fit model are shown in Fig.|2] 
and the line fluxes are listed in Table [l] The uncertainty in the 
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Fig. 3. Kinetic temperature profile for the circumstellar enve- 
lope of W Hya derived from our best-fit CO model. 

Table 1. Modelled line fluxes for CO compared with the obser- 
vations (the unit is K km s _1 , except for the ISO lines where the 
unit is WcrrT 2 ). The observed data are taken from the follow- 
ing references : a - Kerschbaum & Olofsson ( 1999), b - JCMT 
archive, c - Cernicharo et al. ( 1997 1, d - Olofsson, private com- 
munication, e - Young ( 1995 1, f - Barlow et al. 1 1996 1 



Transition 


Instr. 


Observed flux 


Model flux 


Ref 


1-0 


SEST 


0.82 


0.66 


a 


2-1 


JCMT 


18.7 


11.8 


b 


2-1 


IRAM 


22.0 


26.2 


c 


3-2 


JCMT 


40.2 


39.8 


d 


3-2 


CSO 


29.0 


22.7 


e 


4-3 


JCMT 


44.6 


60.8 


b 


4-3 


CSO 


40.0 


34.8 


e 


16-15 


LWS 


3.3 10- 20 


2.4 10~ 20 


f 


17-16 


LWS 


2.3 10- 20 


2.4 10- 20 


f 



mass-loss rate from fitting the line profile is of the order of 
50% within the adopted circumstellar model (see Schoier & 



Olofsson 2001 ). The final model also gives the LSR velocity 
for W Hya as 40.7 kms -1 , and an envelope gas expansion ve- 
locity of 7.0 km s _1 . These values are consistent with those de- 
rived from the Odin H2O data. 

The final kinetic temperature profile is shown in Fig. [3] We 
note that our derived temperature is between 100-150 K higher 



than that obtained by Zubko & Elitzur (2000 1 in the region 
where the H2O molecules are located (r < 2 10 15 cm). This 
is not inconsistent considering our lower mass-loss rate. 

3.2. H 2 lines 

We model the water lines observed by Odin, ISO-SWS 
(Neufeld et al. fT996l . and ISO-LWS (Barlow et al. 1996 1. The 
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code used for this differs from the code used for CO. The main 
reason for a different treatment is the extremely high optical 
depths in the water lines preventing the convergence of the 
Monte-Carlo calculation. The technique employed here is the 
accelerated lambda iteration (Rybicki & Hummer 1991 1, which 
is used to calculate the excitation in the 45 lowest rotational 
levels of ortho- and para-HaO. 

We use the collisional rates calculated for H2O and He 
(Green et al. 1993 ), with a scaling factor of 1 .4 to correct for the 
mass difference between H2 and He. However, a more recent 
calculation for the collision between H2O and ortho- and para- 
H2 by Phillips et al. (1996) suggests that adopting the colli- 
sional rate with He and using the simple scaling factor can give 
very different rates. Unfortunately, the collisional rates with H2 
were calculated up to 140 K, while the temperature in the cir- 
cumstellar envelope of an AGB star can be up to 1 000 K. If the 
H20-He rate is underestimated by a large factor, it may have a 
repercussion on the calculated line fluxes as the water lines are 
sub-thermally excited in the envelope of W Hya. 

The envelope gas expansion velocity, and the kinetic tem- 
perature law used are those obtained from the CO line mod- 
elling. The outer radius of the water molecules is taken from 
the measurements of the OH masers by Szymczak et al. ( 1998 1, 
who give an outer OH shell radius of 1.3", i.e., 1.5 10 15 cm at 
the assumed distance. 

The mass-loss rate and the ortho-H20 abundance are free 
parameters in the fitting procedure where the integrated line 
intensities of the Odin and ISO lines are used as observa- 
tional constraints to find the best-fit model in the% 2 sense. The 
observed integrated line intensities and their uncertainties are 
given in Table|2] We obtained ISO fluxes by fitting a Gaussian 
to archived data with the latest calibration file (OLP version 
10), rather than using the line fluxes from Barlow et al. i 19961 1 
and Neufeld et al. ( 1996 1. The error in the observed flux re- 
flects only the estimation of the baseline and does not include 
possible systematic error or that of the absolute flux. 

The results of the line emission modelling are given in 
Fig. 0] in terms of the minimum x 1 f° r a given combination 
of ortho-H20 abundance and mass-loss rate. Clearly, the cir- 
cumstellar H2O line emission suggests a relatively low mass- 
loss rate. Using a l<x uncertainty we infer a mass-loss rate 
of (2.5±0.5) 10" 7 M o yr _1 , and an ortho-H 2 abundance of 
(1.0+0.4) 10 3 . The best-fit model line profile for the Odin 
557 GHz line is presented in Fig. [5] superimposed on the ob- 
served spectrum. The resulting Odin and ISO lines are shown 
in Fig. [6] although at a much higher spectral resolution than 
achieved by ISO. 

We have also modelled the para-water lines present in the 
ISO spectrum. Since the ortho- and para-water molecules do 
not interact with each other, they can be modelled indepen- 
dently. In order to fit these lines for the mass-loss rate derived 
above, the required para-water abundance is 1.1 10~ 3 (with an 
estimated uncertainty similar to that of the ortho-H20 abun- 
dance), The ortho-to-para ratio is close to unity, in agree- 
ment with that reported by Barlow et al. ( 1996). This ratio is 
expected to be three when considering newly formed H2 on 
the grain surface at high temperature (Tielens & Allamandola 
1987 ). The ratio drops to zero for a cold gas. 
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Fig. 4. The minimum x 1 f° r a given combination of ortho- 
H2O abundance and stellar mass-loss rate (the number at each 
point). The Odin and ISO LWS/SWS water line fluxes have 
been used as observational constraints (see Table|2j(. 



0.06 - 



0.04 - 



0.02 - 




-0.02 



v lsr (km/s) 

Fig. 5. A comparison between the Odin H2O 1 10 - lot observa- 
tions (histogram) and the model profile from the best-fit model. 
The model profile shows clearly the self-absorbed blue wing. 



The total (ortho + para) water abundance according to this 
analysis, consistent with the observed intensities, is therefore 
(2.0+1.0) 10~ 3 . 

The marked increase in the y 1 value for higher mass-loss 
rates is primarily an effect of the increasing difficulty to fit the 
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Table 2. Resulting ortho-f^O line fluxes for the best-fit model 
with a stellar mass-loss rate of 2 10~ 7 Moyr _I and an ortho- 
H2O abundance of 1.1 10~ 3 . The lines listed are the Odin 
line (538.29 ^m) and some ISO lines previously detected by 
a Barlow et al. (I1996> and b Neufeld et al. i l 19961 1. See text for 
details on the observed line fluxes and their uncertainties. 



A (pm) 


Instr. 


Line flux (10 


- 20 WcirT 2 ) 






Obs 


Mod 


538.29 


Odin 


0.36+0.1 


0.45 


179.53" 


LWS 


10.3+1.7 


11.6 


174.63° 


LWS 


11.9+1.2 


8.33 


108.07" 


LWS 


22.8+4.0 


15.3 


66.44" 


LWS 


32.7+12.1 


21.5 


40.69* 


SWS 


23.4+5.3 


32.6 


31.77* 


SWS 


36.3+3.7 


31.1 



SWS lines. The model line fluxes are significantly lower than 
the observed ones due to the fact that the blue edge of the pro- 
file is in absorption, cancelling out the emission from the red 
side of the profile (see Fig |6}. The absorption is not so pro- 
nounced in the low-mass-loss-rate model, hence the overall line 
fluxes agree more with the observations. It is unfortunate that 
the resolution of the ISO lines is not high enough for the lines 
to be resolved. It should be remarked here that the SWS lines 
come from the very inner region of the circumstellar envelope, 
and our modelling of these lines is more uncertain. Fitting only 
the LWS lines gives a weaker constraint on the combination of 
the mass-loss rate and water abundance. 

We have also investigated the effect of the assumed ra- 
dius of the water envelope by varying this with respect to that 
adopted from the OH maser data. By increasing the radius by 
a factor of three, the resulting line fluxes changed by less than 
50% for our best-fit model. The line fluxes are primarily ex- 
citation limited and they depend only weakly on the adopted 
H2O envelope size. Netzer & Knapp I1987> derived a size of 
the water envelope from a theoretical model of H2O photodis- 
sociation 

R oul = 1 10 15 (M/ICT 7 ) - 7 cm, (2) 

where the mass-loss rate is given in M yr~'. We find that the 
radius obtained from the OH maser data is fully consistent with 
this relation for the mass-loss rate of 2.5 10~ 7 M yr~' obtained 
from the H2O line flux fitting. 

In our calculation, we did not include the effect of radiative 
excitation through the vibrational lines. W Hya has a deep H2O 
absorption band at 3 /vm, due to the stretching mode, as seen 
in the ISO SWS spectrum (Justtanont et al. 2004 1. This band 
occurs close to the stellar photosphere (r ex ~ 10 3 K) and most 
likely a significant fraction of the stellar flux at the wavelengths 
of the vibrational lines is absorbed in it (Ryde et al. 2002). The 
detection of H2O maser emission from rotational lines in the vi 
vibrational state towards W Hya shows that radiative excitation 
to this state occurs (Menten & MelnickQ989|l. Deguchi & Rieu 
( 1990 1 found that vibrational radiative excitation plays a minor 
role when calculating the rotational line fluxes in the ground 
state. On the contrary, Troung-Bach et al. ( 1999 1, in their study 



of the O-rich Mira variable R Cas, found that the ground-state 
line fluxes decreased by about a factor of two when radiative 
excitation through the V2 vibrational state was removed. We be- 
lieve that the uncertainty in the H2O abundance introduced by 
omitting radiative excitation through vibrational states is less 
than our estimated uncertainty from the line modelling. The ef- 
fect is likely to lower the abundance somewhat. We note here 
that Truong-Bach et al. derive an H2O abundance for R Cas 
which is remarkably low (in terms of our modelling results), 

1 10~ 5 , considering that R Cas and W Hya have very similar 
mass-loss rates, luminosity, distances, and ISO line fluxes. 

4. The mass-loss rate of W Hya 

There is a large discrepancy when it comes to the determina- 
tion of the mass-loss rate of W Hya as a wide range (two or- 
ders of magnitudes) of the mass-loss rate has been reported 
over the years. In this paper we have shown that the circum- 
stellar CO radio lines suggest a low mass-loss rate, about 

2 lO _7 M yr _I . This result has a rather low uncertainty pro- 
vided that the adopted circumstellar model and CO abundance 
(2 10~ 4 ) are correct. Our result is also consistent with those of 
Wannier & Sahai ( fT986l . Young (IT9951 . and Olofsson et al. 
(2002). Once these values are scaled to the same distance, ve- 
locity, and CO abundance, the range of the mass-loss rates is 
(S^lO^Moyr- 1 . 

Previous reports on the mass-loss rates derived from H2O 
line emission are based on ISO results. The mass-loss rate esti- 
mates range from as high as 10~ 5 (assumed distance 95 pc and 
SWS lines; Neufeld et al. fT996t to 6 10" 7 M Q yr 1 ( assumed 
distance 130pc and LWS lines; Barlow et al. 1996). Barlow 
et al. used a water abundance of 8 10~ 4 for r < 6 10 14 cm and 

3 10~ 4 beyond this radius. Neufeld et al. do not give their as- 
sumed water abundance, but it must have been low, according 
to our results (see Fig.0). We find that the Odin and ISO H2O 
line intensities suggest a mass-loss rate of 2.5 10~ 7 M yr~', a 
value fully consistent with our CO line modelling result. 

The OH maser line data give some hints on the mass-loss 
rate. Recently, the 1612 MHz line was detected (Etoka et al. 
2003 ), and it is substantially weaker than the main lines at 1665 
and 1667 MHz. This is an indication of a low mass-loss rate 
(Szymczak & Engels 1998). Likewise, as pointed out above, 
the measured size of the OH maser envelope is consistent with 
a low mass-loss rate. 

Justtanont et al. ( 2004) derived a low dust mass-loss rate by 
fitting the SED. Using their derived dust mass-loss rate com- 
bined with a dynamical calculation of a dust-driven wind, we 
infer a total mass-loss rate of 1.4 10~ 7 M yr _1 . This value is 
only a lower limit because the different dust components used 
in the SED modelling do not account for the 13/vm dust fea- 
ture. Once the carrier is found for this feature and taken into 
account in the model, it is possible that the dust mass loss rate 
can increase by up to a factor of two. The indication from the 
dynamical calculation is W Hya has a relatively low mass-loss 
rate. 

However, Zubko & Elitzur (2000) reported a rather high 
mass-loss rate when combining SED modelling and a dynam- 
ical model for the outflow, 2.3 lO~ 6 M yr _1 (at a distance of 
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1 15 pc). The high mass-loss rate is because the authors did not 
fit the 10/vm silicate feature which is the key in deriving the 
dust mass loss rate. 

Based on these results we conclude that there is strong evi- 
dence that W Hya has a rather low mass-loss rate, of the order 
a few lO _7 M yr _1 . 

5. The circumstellar H 2 abundance 

In this work, we derive a consistent mass-loss rate from both 
the CO and H 2 modelling of (2.5+0.5) 10 7 Movr" 1 . The 
observed line fluxes are reasonably well produced. Our de- 
rived water abundance for this mass-loss rate is exceptionally 
high, (2.0+1.0) 10 -3 . Additional uncertainties in this estimate 
are provided by the effects of radiative excitation through vibra- 
tional states, and the uncertain values for the collisional cross 
sections. 

This derived total water abundance violates the limit set 
by the cosmic abundance of atomic oxygen. In the stellar 
equilibrium chemistry of an oxygen-rich environment, all the 
available carbon should be locked up in the form of CO, 
and the rest of the oxygen goes into the formation of H2O 
molecules, with a very small fraction bound into grains and 
other molecules. The cosmic abundance of Carbon and Oxygen 
relative to H is 3.6 10~ 4 and 8.5 10~ 4 , respectively (Anders & 
Grevesse 1989). Therefore, unless there is a separate influx of 
water from another source the water abundance relative to H2, 
/ H ,o = [H 2 0]/[H 2 ], cannot exceed 10~ 3 . 

The extra water mass required to explain the abundance 
needed to model the observed water line emission is 

Mh 2 o = M(/ H ,o - 10~ 3 )(rH,o/v)(mH 2 o/mH 2 ) ( 3 ) 



where n-j,o is the radius of the water shell, and v the gas expan- 
sion velocity. The total extra mass needed is (6.1+3.0) 10 26 g, 
i.e., (0.1±0.05)M e . This is considerably smaller than the esti- 
mated mass of water ice in the Solar System, assuming most 
of the present mass of water is bound in comets. The current 
estimate of the original mass of the Oort Cloud is highly uncer- 
tain and varies between 14-1 000 M e (Bailey & Stagg 1988; 
Weissman ll99U Marochnik et al. 1988). We speculate that icy 
bodies can be either part of the disk (and planets) left over af- 
ter the formation of the central star, or a remnant envelope of 
the original cloud from which the star formed, reminiscent of 
the Oort Cloud. The disk/envelope remained in stable orbits 
until the star ascended the AGB. Because of the much higher 
luminosity of the AGB star, compared to that in its main se- 
quence phase, outgassing from icy bodies will contribute to 
the observed high water abundance. Of course, this requires a 
steady supply of new water, since photodissociation effectively 
destroys the water molecules on a short time scale. A simi- 
lar hypothesis has been put forward to explain the presence of 
gaseous water around the carbon star IRC+10216 (Melnick et 
al. l200Tl . 



6. Summary 

Observations of the ground-state transition of ortho-water, 1 10- 
loi, in the circumstellar envelope of W Hya were done us- 
ing the Odin satellite. The line is centered at Vlsr = 42kms~ 1 
and its width corresponds to a gas expansion velocity of about 
7kms _1 . These results compare well with those of previous 
CO observations (e.g., Kerschbaum & Olofsson 1999), as well 
as the new CO data presented in this paper. The line flux is 



K. Justtanont et al.: W Hya through the eye of Odin 



7 



(3.6+1.2) 10 21 W crrT 2 , consistent with the SWAS observa- 
tion of the same line (Harwit & Bergin 2002}. 

To fit the Odin and ISO water lines, the ortho-water abun- 
dance and mass-loss rate are required to be (1.0±0.4) 10 -3 
and (2.5+1.0) 10 -7 M yr _1 , respectively. Similarly, the para- 
water abundance is estimated to be (1.1 ±0.8) 10 3 by fitting 
ISO lines. This mass-loss rate is consistent with that derived 
from our CO line modelling, as well as that obtained by com- 
bining SED modelling and a dynamical model of the outflow. 
From our calculation, we show that the line shapes of these wa- 
ter transitions are different. Unfortunately, from ISO, we have 
only the information on the line intensities. It is noted here that 
the Heterodyne Instrument for the Far-Infrared (HIFI) aboard 
the Herschel Observatory which will be launched around 2008 
will have a much higher resolution which can resolve some of 
these ISO lines, thereby giving another constraint on the model 
parameters. 

If the mass-loss rate for W Hya is as low as (2.5+0.5) 10~ 7 
M Q yr _1 , the total water abundance needed to explain the ob- 
served line fluxes is (2.0+1.0) 10~ 3 . However, if we assume the 
cosmic abundance for carbon and oxygen, the maximum total 
water abundance, relative to H2, cannot be higher than 10 -3 . 
A way to reconcile the high total water abundance needed to 
model the H2O lines of W Hya is an extra injection of water 
in the system, possibly from evaporation of icy bodies, such 
as planets or comets. This is the suggested explanation for the 
detection of water in IRC+10216, the IR-brightest carbon-rich 
AGB star (Melnick et al. 2001 ). The excess H2O mass needed 
to explain the Odin and ISO observations in W Hya is small, 
of the order of (0. 1 ±0.05) M ffl , compared to the estimate of the 
original mass of the Oort Cloud. 
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